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Abstract. The observation of the motion of particles and light near a gravitating object is until now the only way to explore 
and to measure the gravitational field. In the case of exact black hole solutions of the Einstein equations the gravitational 
field is characterized by a small number of parameters which can be read off from the observables related to the orbits of test 
particles and light rays. Here we review the state of the art of analytical solutions of geodesic equations in various space-times. 
In particular we consider the four dimensional black hole space-times of Plebahski-Demiahski type as far as the geodesic 
equation separates, as well as solutions in higher dimensions, and also solutions with cosmic strings. The mathematical tools 
used are elliptic and hyperelliptic functions. We present a list of analytic solutions which can be found in the literature. 
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INTRODUCTION 

In the context of General Relativity the gravitational field is the curvature of space and time. This curvature cannot be 
measured directly but has to be inferred from observations. As gravitational waves are not directly detected so far, the 
only way to measure the gravitational field is by the observation of the motion of massive objects and light moving 
in the gravitational field Most suitable for this task are test objects including (i) neutral point particles which are 
described by a normalized four velocity vector which obeys the standard geodesic equation, (ii) light rays, which are 
null vectors also obeying the geodesic equation, (iii) light including polarization, (iv) particles with spin obeying the 
Mathisson-Papapetrou equations 12, Ht], (v) extended particles with spin and mass multipoles which are described 
by the Mathisson-Papapetrou-Dixon equations |01, (vi) electromagnetic fields which are described by the Maxwell 
equations minimally coupled to the gravitational field, and (vii) scalar and spinorial fields obeying the minimally 
coupled Klein-Gordon or Dirac equation. For many applications in astronomy it is sufficient to consider the motion of 
point like particles and light rays. In some cases it is necessary to take the spin and the mass multipoles of test bodies 
into account. Beyond the test object approach self force and radiation reaction effects have to be considered js!]. Here 
we are considering the motion of point like test bodies and light rays only. 

As the motion of test particles and light is of great importance it is most desirable to derive analytic expressions for 
geodesics. Beside the scientific value of analytical solutions itself, first of all analytic expressions enable a systematic 
study of all possible solutions and of the complete parameter space. The structure and characteristics of the solutions 
can be explored. This includes the derivation of observable effects like the perihelion shift or the light deflection. 
Furthermore, analytical solutions are a starting point for approximations in complex scenarios where analytical 
methods fail or for post-Newtonian approximation. They may also serve as test beds for numerical methods. 

In this contribution we review analytical solution methods for the geodesic equation in various space-times beyond 
solutions in terms of elementary functions. We describe how solutions in terms of elliptic or hyperelliptic functions 
can be derived and provide a guide to the literature, which addresses such problems. 


GEODESIC MOTION 


In General Relativity the motion of a point like neutral test particle is described by the geodesic equation 


ds^ 


+r^ 


dx^ dxP _ 
ds ds 


r^vp — 2 ^^°^ if^vgpa + dpgva — dagvp) 


( 1 ) 


where 


( 2 ) 



are the Christoffel symbols, which depend on the space-time metric {gjiy), and s is an affine parameter. If the metric 
(gjiv) is axially symmetric there are the two constants of motion 

E — (3) 


which are related to the Killing vectors and and can be interpreted as the energy and the specific angular 
momentum in direction of the symmetry axes. Here u denotes the four-velocity. In addition, we may normalize 
gfiydx^dx'' = £, where e = 0 for light and e = 1 for massive particles. However, these three constants are in general not 
sufficient to separate the geodesic equation O- If we have even spherical symmetry we may restrict the motion to the 
equatorial plane but this is not possible in general. Surprisingly, a fourth constant of motion, which is not connected to 
an obvious symmetry, can be found for a quite wide range of space-times. In 1968 Carter (@] found such an additional 
constant of motion, which can be derived as a separation constant. As a result, the geodesic equation in a wide range 
of space-times can be rewritten in the form 


dx\^ 
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where x is (a function of) the radius r or latitude 0 , P{x) is a polynomial, and X is an affine parameter, which in general 
does not coincide with s. (Note that in spherically symmetric space-times often A can be replaced by the angle (p.) The 
equations for the angle (p and the coordinate time t can then be derived from Q and take the form 

^=f{r)+g{e), (5) 


where xis (p or t and /, g are rational functions. 

For example, timelike geodesic motion and even the motion of charged particles in Kerr-Newman-Taub-NUT-de 
Sitter space-time 10,01 may be reduced to the form (IHi, (|5]|. For lightlike motion in these space-times in addition 
accelerated gravitating objects may be considered 10 ]. Also, timelike and lightlike geodesic motion in some higher 
dimensional generalizations of Kerr-Newman-de Sitter space-time may be reduced to the form (IHi, (|5]) as well as Kerr 
and Schwarzschild space-times pierced by cosmic strings and the geodesic motion in some spherically symmetric 
regular black holes. A more detailed account on applications to various space-times can be found in the section on 
analytical solutions in the literature. 


ANALYTICAL SOLUTION METHODS 

The analytical solution methods for differential equations of the type (|4l) depend on the degree of the polynomial 
P. If the degree of P is two or lower than the differential equation can be solved in terms of elementary functions. 
If the degree is three or four the equation can be solved in terms of elliptic functions. If it is of even higher degree 
hyperelliptic functions are used to solve the equation. The solution methods in terms of elliptic and hyperelliptic 
functions are outlined in the section, where we restrict ourselves to degrees below six and only shortly mention the 
even more general case. 


Elliptic functions 


Differential equations of the type 

^ =P3Aix), Jc(yo)=.*o (6) 

where F 3 4 is a polynomial of degree three or four, can be solved in terms of elliptic functions. These are meromorphic 
functions with two linearly independent periods 2c0i , 20)2 £ C, i.e. ^ ©2 ^ A non constant elliptic function has a 
finite number of poles in each period parallelogram {z£C|30<f< 1: z = zo+2tcoi +2t(02}, whose sum of residues 
has to vanish (see e.g. SIH). Historically two different approaches to elliptic function were developed named after 
Jacobi and Weierstrass. 



The approach by Weierstrass is based on the most simple elliptic function which has a double pole in zero with 
vanishing residue. This function, called the Weierstrass p function, can be constructed as a series 
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where the last term has to be added to achieve absolute and uniform convergence. By counting poles and using that 
elliptic functions without poles are constants it can be shown easily that the Weierstrass p function solves the initial 
value problem (see e.g. 


^ =Ax^-g2x-g^^ x( 0 )=oo ( 8 ) 

where the invariants g 2 , g3 are given by g 2 = 60L(„_m)^(op)(2na)i +2ma)2)^^ g3 = 140L(„,m)^(o.O)(2n®i + 2 ma) 2 )^®. 
Note that they are connected to the zeros ej, J = 1,... ,3 of Pfv(x) := 4x^ — g 2 X — g 3 by 

ei +62 + 63=0, -4(6162 + 6163+ 6263 ) =^ 2 , 4616263 = g 3 . (9) 

For the applications we have in mind here usually fV and, therefore, the zeros 6 / and the invariants g 2 , g 3 are known 
instead of the periods 2a)i,2- The periods can be calculated from the invariants by integration: in general, for two 
different closed integration paths 71 2 , 
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as the zeros Cj are singularities of the integrand. All periods (2na)i +2m(02), n,m € Z, can then be determined by 
integrating along all possible paths. The basic periods 2coi 2 can e.g. be determined by 
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if all 6; are real and ordered as 61 < 62 < 63 or if 61 =62 and 63 G M, where the paths are assumed not to cross each 
other. (Here we used the convention that (£)\ G K; other conventions use (£)\ = ©2 if 61 —62 ) 

The series expansion (|2l) is not very well suited to actually compute the Weierstrass p function as it converges quite 
slowly. Alternatively it can be calculated iteratively ifTTI] . using a power series with Eisenstein series as coefficients, or 
in terms of theta functions l&lli. The Riemann theta function is defined as 


0 [tv + w](z; t) = exp(;r/(m +v)'(T(m + v) +2z + 2w)), 
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where z G C'^, T is a g x g symmetric matrix with positive definite imaginary part. Tv + w G is the characteristic, 
and g is the genus, here g = 1 . 

The general differential (| 6 ]l can be solved in terms of the Weierstrass elliptic function by converting the problem to 
the standard form (l 8 ]l. If F’ 3,4 is of degree four first apply the substitutions: = +xp where P 3 ^ 4 {xp) = 0 to convert to 
a polynomial of degree three. Subsequently, or if P 3 4 was of degree three in the first place, substitute ^ = ^(4z+ y) 
if P = Y.it'ix' is the polynomial of degree three. Then 
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The approach by Jacobi is based on the inversion of elliptic integrals. The sine amplitude function sn can be defined 


as 
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where 0 < ^ < 1 is the modulus, w G [0,1], and z € R. It has poles in iK{l<') and 2K{k) + iK{k') where K{k) =F{l-,k) 
is the complete elliptic integral of the first kind and {k')^ = l —k^. The sine amplitude sn can also be written in terms 
of the Riemann theta function. Indeed, it is often defined in this way. 


sn(z;k) = k 2 
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where T = and u = 

The general differential equation (| 6 l) can be solved in terms of Jacobian elliptic functions by applying a substitution 
which converts the problem to the form 




(16) 


This substitution depends on the degree of Rs 4 and its number of complex zeros as well as on the type of orbit you 
want to obtain. For a list see e.g. 11 121] . 


Solutions in terms of hyperelliptic functions 

In the previous section we mentioned that both the Weierstrass and the Jacobi elliptic functions can be defined in 
terms of theta functions. To generalize the results on elliptic functions to differential equations of the form 
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where 7*5 g is of degree five or six, we first need to define more general functions. In (fTOl l we found that all periods of 
the solution function of (l 6 ]l can be determined by different integration paths. As there the integrand had four singular 
points we had two independent periods. Here, the corresponding integrand has six singular points and, therefore, we get 
four independent periods along four integration paths ai 2 , ^ 1 . 2 , which may be chosen such that integrals over a\ 2 are 
real. However, a function on C cannot have more than two independent periods. This means that we have to consider 
more general functions. Let us first introduce the holomorphic differentials dzi and the meromorphic differential 
by 
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where P{t) = 4 t^s+i is a polynomial of degree 2g + 1 and g is the genus. We may then define the 

corresponding g x g period matrices of the first kind (2a)i ), { 2 ( 02 ) and of the second kind (2rji), ( 2772 ) by 
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Here ai,...,ag, b\,...,bg are the 2g independent closed integration paths. If P has only real zeros ej ordered as 
Cj < Cj+x one may e.g. choose aj as the path from e 2 j-i along the real axes to e 2 j and back (with reversed sign of 









the square root) and the path bj from eij along the real axes to eij+i and back (again with reversed sign of the square 
root). 

To generalize the solution methods outlined in the previous section we first need to consider the Jacobi inversion 
problem 
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Note that forg = 1 and 02 = 0 we recover®. The g solutions .ry(yi,... of (ISlT i can be given in terms of generalized 
Weierstrass functions. These are defined by the theta function via the Kleinian sigma function 
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where z € C^, ic — ri(2c0i) \ and C is a constant with does not matter in the present context. The vector of Riemann 
constants with base point at infinity, Kx, S C'^ can be calculated by 
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where e, are the zeros of the polynomial P. This means that is a sum of normalized half periods, in particular 
Koo = t(|, j)' + (0, 4)^ if g = 2. Note that T = a)j“* 0)2 is always symmetric with positive definite imaginary part. For 
further details see e.g. lll3ll . The solutions of (ISTTi are then given by the solutions of 

+ Y,Pgiiyi,---:yg)x‘^^ =0- ( 25 ) 

(=1 

Let us turn back to the differential equation (fTTl i and genus g = 2. If is of degree six the equation can be 
reformulated as = ^/P${^ by a substitution x = +xp, where xp is a zero of P^ e P^ is a polynomial of 
degree five. Such a differential equation, or, if P 5 g was of degree five in the first place, can then be rewritten in the 
form 


f' = t{yo)=to, (26) 

dy 

where i = I or / = 2, by an appropriate normalization. We may then consider this equation as part of the Jacobi 
inversion problem. 
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We may then find the solution to (l26T l as the limiting case X 2 ^°° of the Jacobi inversion problem (|2TI) in the following 
way: first observe that 
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where Oj(z) is the derivative of o with respect to Zi- From the comparison of d^ . where g = 2 and X 2 with (l27l) 
we may identify y, in (l28l) with our physical coordinate y, in (l27l) . Fortunately, we get automatically rid of the other yj, 
j ^ i, by the same limiting process X 2 —?■ 0 °. This is because the set of zeros of the theta function z 0[AL]((2a))^^z), 
which is a one dimensional submanifold of C^, is given by all vectors z = (zi,Z 2 ) which can be written as z; = f^dzi 










with the same x and dzi as in (fTsT l. see e.g. iB. This is exactly true for the vector {yi,y 2 ), which means that we may 
write yj = f{yi) for some function /. As the zeros of the theta function are also zeros of a we can simplify (l28T l to 

t = -^{yuf{yi)) or t = -^{f{y 2 ),y 2 )- (29) 

CT2 <72 

with y, as defined in (l27l i. 

Note that the solution method for differential equations of the type (fTTI) can be generalized even further. As long as 
we have on the right hand side of (fTTl i a polynomial of any degree, we may consider the equation as part of Jacobi’s 
inversion problem. As an additional complication we have to find g — 1 condition on the y, ’s on left hand side in (|2T]) 
to restrict the problem to a one dimensional submanifold of C'^. For details see B. 

ANALYTICAL SOLUTIONS OF GEODESIC EQUATION IN VARIOUS SPACE-TIMES 

In this section we will collect applications of the methods outlined in the section on geodesic motion to various 
space-times. For older literature we refer to Sharp ifl^ who collected most of the papers on geodesic motion in Kerr- 
Newman space-time and subclasses, which were available at that time. Partly this is still quite complete but we also 
try to update his collection (with respect to analytical solutions). Note that we only consider analytical solutions to 
general timelike and lightlike geodesics (with an electric or magnetic charge, as applicable). In particular, we do not 
list the vast literature on equatorial motion in Kerr space-time. Of course, we do not claim that our list is complete. 


Particle motion in standard black hole space-times 

Schwarzschild 

Regarding analytical solution methods the list of Sharp already contained the complete set of solutions. Most 
notably, this includes the works by Hagihara 11^ . who derived the analytical solutions in terms of Weierstrass elliptic 
functions, and Darwin iflslfT^ . who used Jacobian elliptic functions and integrals. 

Accordingly, observables like the perihelion shift, the scattering angle of flyby orbits, the light deflection (which is 
the flyby orbit of light) can be described in terms of (complete) elliptic integrals. 


Reissner-Nordstrom 

Surprisingly, the analytical solutions to the geodesic equation in Reissner-Nordstrom space-time seem to be consid¬ 
ered first only in 1983 by Gackstatter although it can be handled completely analogously to the Schwarzschild 
case. He studied bound timelike geodesics and light in terms of Jacobian elliptic integrals and functions. Recently, 
Slezakova gave a comprehensive analysis of arbitrary timelike, lightlike, and even spacelike geodesics. Grunau 
and Kagramanova ll^ solved the equations of motion of electrically and magnetically charged particles in Reissner- 
Nordstrom space-time in terms of Weierstrass elliptic functions. Gackstatter and Grunau/Kagramanova also considered 
the perihelion shift and the deflection of light, which can be given in terms of (complete) elliptic integrals. 


Taub-NUT 

Timelike geodesics were studied by Kagramanova et al ll^ in terms of Weierstrass elliptic functions. One particular 
feamre of the geodesics in NUT space-times is that the orbits do not lie on an orbital plane but, instead, on an orbital 
cone with an opening angle which depends on the NUT parameter, which can also be considered as a gravito mag netic 
mass. In particular, the issue of geodesic incompleteness at the horizons as discussed by Misner and Taub 11^ . and 
the analytic extension of Miller, Kruskal and Godfrey ll^ was addressed. 


Kerr 


Most of the older literature on Kerr space-time is concerned with the much simpler particular case of equatorial 
geodesics. We refer to Sharp ifT^ here for these early works. Note that in terms of the proper time (or the corresponding 
affine parameter for light) the equations of motion are still coupled. Therefore, most of the analytical solutions before 
the introduction of the Mino time implicitly included integrals over the latitude or the radius, see e.g. Kraniotis 
llz^ or Slezakova ll^ for a review. As notable exception, Cadez et al ll2^ introduced already in 1998 a similar 
parameter (called P, see their equation (34)) as they considered the motion of light. 

After the introduction of the Mino time, in 2009 Fujita and Hikida ll2^ used this new affine parameter to derive 
the analytical solution for bound timelike geodesics in terms of Jacobian elliptic functions. They also elaborated the 
expressions for perihelion shift and Lense-Thirring effect in terms of Jacobian elliptic integrals. General timelike 
geodesics and lightlike motion were treated shortly after tha t by Hackmann in 2010. Recently, also an analytical 
expression for the gravitomagnetic clock effect was derived IBHI . 

Note that Kraniotis ll^ also derived analytical solutions for lightlike geodesics in terms of hypergeometric func¬ 
tions. 


Kerr-Newman 

Charged particle motion was considered by Hackmann and Xu in terms of Weierstrass elliptic functions. Note 
that also a magnetic charge of the black hole was included. Due to many parameters involved the number of orbit 
configurations is very large. This includes orbits crossing the horizons or r = 0. For the sake of completeness, also 
black holes endowed with magnetic charge has been considered. This has a big impact on the orbital types: Not only 
does the motion deviate from the symmetry to the equatorial plane but also stable off-equatorial circular orbits outside 
the horizon do exist in this case, which are not possible elsewhere. 


Schwarzschild-de Sitter 

This space-time is also named after Kottler. Note that on the level of the differential equation, lightlike geodesics 
in Schwarzschild-de Sitter space-time are identical with the lightlike equations of motion for Schwarzschild, as the 
cosmological constant can be absorbed in the definition of just a single parameter. Analytical solution are given e.g. in 
Gibbons et al ll^ . General timelike geodesics in Kottler space-time can be treated in terms of hyperelliptic functions 
as elaborated by Hackmann and Lammerzahl ll^ . 

The motivation for searching for an analytical solution of the geodesic equation in Schwarzschild-de Sitter space- 
time was to investigate the suggestion being made that the Pioneer anomaly could originate from the influence of the 
cosmological constant on the motion of test bodies. With the above analytic solution it could be shown that the effect of 
a cosmological constant is of the order of cm only and, thus, cannot be responsible for the the anomalous acceleration 
which has been unambiguously measured. In fact, later the anisotropic thermal radiation could be identified as origin 
of the anomalous acceleration ifj^ . 


Reissner-Nordstrom-de Sitter 

The equations of motion for general timelike geodesics were solved in ll^ . The motion of photons was very 
recently analytically calculated by Villanueva et al 13^ for a negative cosmological constants using Weierstrass elliptic 
functions. 


Taub-NUT-de Sitter 

The transition from Taub-NUT to Taub-NUT-de Sitter is similar to the transition of Schwarzschild to Schwarzschild- 
de Sitter. For a not too big cosmological constant an additional cosmological horizon appears which enlerges the 
possible orbital types in an obvious way. The timelike motion in this case has been analyzed in 


Kerr-de Sitter 


The equations of motions for timelike geodesics were analytically solved by Hackmann et al MM in terms of 
hyperelliptic functions. Note that Kraniotis S also derived analytical solutions for lightlike geodesics in terms of 
hypergeometric functions. 


Kerr-Newman-Taub-NUT-de Sitter 


The general solution for timelike geodesic was shortly outlined in 14^ . Analytical expressions for the perihelion 
shift and the Lense-Thirring effect were derived in ll43l . 


Particle motion in higher dimensional space-times 

Spherical symmetry 

Timelike and lightlike motion in higher dimensional Schwarzschild, Reissner-Nordstrom, Schwarzschild-de Sitter, 
and Reissner-Nordstrom-de Sitter was considered in ll^ and more generally in ifisll . 


Myers-Perry 

Myers-Perry space-times are higher dimensional generalizations of the 4-dimensional Kerr space-time. Kagra- 
manova and Reimers ll4^ derived analytical solution to the geodesic equation in a five dimensional Myers-Perry 
space-time assuming that the two parameters of rotation coincide. Their results are valid for both massive test particles 
and light. Also the various observables, that is, the two types of perihelion shift as well as two types of Lense-Thirring 
effect have been derived. 


Black rings 

In higher dimensions black holes may have a more complicated horizon topology than in four dimensions. In 
five dimensions, the horizon may have the topology of a ring which are accordingly called black ring space-times. 
Grunau et al considered geodesics in a singly spinning black ring spacetime in five dimensions. In general the 
Hamilton-Jacobi equation is not separable for these space-times. For the special cases where it is separable Grunau et 
al derived analytical solutions for timelike and lightlike motion. The same holds for the doubly spinning and charged 
five dimensional black ring considered by Grunau et al in 2013 ll4^ . 


Black strings 

Grunau and Khamesra considered particle motion in (rotating) black string space-times, which are given as 
generalizations of Schwarzschild and Kerr space-times including a compact fifth dimension. They derived analytical 
solutions for timelike and lightlike motion in both the non-rotating and the rotating space-time. 


Particle motion in cosmic string space-times 

Schwarzschild black hole with cosmic string 

Owing to the possible connection to string theory, cosmic strings have gained a lot of renewed interest over the 
past years. They are topological defects that could have formed in one of the numerous phase transitions in the early 


universe. Therefore it is of principle interest to consider the possibility of black holes pierced by an infinitely thin 
cosmic string and to work out possible observational consequences. 

This has been catTied through in lli^ . In this paper the geodesic motion in the space-time of a Schwarzschild black 
hole pierced by an infinitely thin cosmic string has been solved analytically for massive bodies and for light rays, 
respectively. The solutions of the geodesic equation depend on the particle’s energy and angular momentum, the ratio 
between the component of the angular momentum aligned with the axis of the string and the total angular momentum, 
the deficit angle of the space-time as well as the mass of the black hole. One main difference to the orbits in the 
Schwarzschild case is that they are not restricted to an orbital plane. The change of the deficit angle affects the radial 
motion significantly. 

In a bound of the energy per unit length of the string was found and compared with experimental tests of 
general relativity. From perihelion shift and light deflection a rough estimate for the string line mass density was 
derived, jx < lO'^^'kg/m. The eventual existence of cosmic strings will also influence the creation of gravitational 
waves since the deficit angle modifies the angular velocity and, thus, the temporal change of the effective quadrupole 
responsible for the emitted gravitational wave. 


Kerr black hole with cosmic string 

For the same reasons there may exist Kerr black holes pierced by a cosmic string. Hackmann et al determined 
the complete set of solutions for this case for massive and mass less particles using the standard techniques developed 
for the Kerr space-time. The solutions of the geodesic equation can be classified according to the particle’s energy and 
angular momentum as well as the mass and the specific angular momentum of the black hole. 

Also the perihelion shift and the Lense-Thiri'ing effect for bound orbits has been discussed. It was proven that 
the presence of a cosmic string enhances both effects JUt]. Comparison of these results with experimental data from 
the LAGEOS satellites yields an upper bound on the energy per unit length of a string piercing the earth which is 
approximately 10*® kg/m. These results also may have relevance to the recently suggested explanation of the alignment 
of the polarization vector of quasars using remnants of cosmic string decay in the form of primordial magnetic held 
loops. 


Particle motion in further space-times 

Thin-shell wormhole 

Diemer and Smolarek recently considered timelike and lightlike motion in traversable Schwarzschild and Kerr 
thin-shell wormholes. 


Regular black holes 

Regular black holes without curvature singularity are considered as an alternative to singular black hole solution 
to avoid the conceptual problems related to the occutTence of a singularity. First regular black holes models were 
constructed by Bardeen in 1968 ifs^ and later provided with a physical reasonable source by Ayon-Beato and Garcia 
Ayon-Beato and Garcia also derived an exact regular solution of the Einstein equations coupled to a nonlinear 
electrodynamics, which may be called the Ayon-Beato Garcia space-time if^ . 

Due to the occurrence of a square root the geodesic equation in this space-time is on the first glance not of the type 
considered here. However, a substitution and the introduction of a new affine parameter analogously to the Mino time 
give a differential equation of the type (fTTI i which can be solved completely if^ . Also the perihelion shift has been 
calculated and its difference from the corresponding effect in a Reissner-Nordstrom space-time has been isolated. 

Though the orbit for a neutral particle could be solved analytically, the orbit for a charged particle leads to a 
mathematically different equation of motion. One obtains 

y'^ + P„{r)y^ + P,„{r) = 0 with y:=f{u)^, 


(30) 


where P„ are polynomials, m is a non-rational function of the radial coordinate (the substitution used before), (p is 
a space-time coordinate, and / is a rational function of u. This equation corresponds to a more general differential 
equation, which may be called a quartic problem according to the corresponding quartic algebraic curve. An analytical 
solution has not yet been found. 


Hofava-Lifshitz gravity 


Hofava-Lifshitz gravity arise from quantum gravity inspired correction to General Relativity. Various black holes 
solutions have been studied where the metric is given by 


ds^ 


N'^{r)dt^ — -^^dr^ — {dd-^ + sin^ ^dtpP') 


(31) 


with = f = I +cir^ — \/cir^ + c^r for some constants ci,C2,C3. In the special case ci = — Avi^, C2 = 0, and 
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C 3 = where Aw is proportional to the negative cosmological constant and a > ^ is an arbitrary parameter 

11 ^ . the geodesic equation can be reduced to the type (|6]l for lightlike motion and to the type (fTTl i for time like motion 
using the substitution u = The analytical solution for this case was derived by Enolskii et al ifs^ . For C 2 7 ^ 0 we 
again encounter a quartic problem. 


OUTLOOK 

In this review we presented the mathematical tools needed to obtain analytical solutions of the geodesic equation 
in a variety of space-times. In a first step, it is possible to solve the geodesic equation in Schwarzschild, Reissner- 
Nordstrom, Taub-NUT space-times with elliptic functions. Owing to the development of the explicit analytic inte¬ 
gration of ordinary differential equations it was possible to solve the geodesic equation also for space-times with 
cosmological constant. This new method can also be applied to further space-times like higher dimensional models, 
space-times related to nonlinear electrodynamics, or to quantum gravity inspired modifications like Hofava-Lifshitz 
gravity. Also particle motion in Gauss-Bonnet gravity models might be treated in this way. This already shows how 
powerful these mathematical tools are. 

However, there are also many other other phenomena and situations where these methods could be applied. As an 
example, the analytical solutions for geodesic motion may be extended to (i) the analytic integration of the Jacobi 
equation lls^ with the application to, e.g., pendulum orbits, cartwheel orbits and helical orbits which may have 
application in satellite geodesy, (ii) the equation of motion for particles with spin in the case that the particle moves 
in the equatorial plane with the spin aligned perpendicular to that plane, (iii) more general for certain motions of 
test particles with mass multipoles. On the level of astrophysical application an analytical timing formula for the 
specific situation of a pulsar orbiting a supermassive black hole may be derived. Here one has to combine the analytic 
description of the motion of a pulsar (treated as test particle) around a massive black hole with the analytic solution for 
the light ray from the pulsar to the observer. This methods may also be applied for the timing between two satellites 
or between a satellite and a clock on Earth. On the long term perspective a further issue is to analytically solve field 
equations in the field of black holes. This has many applications like superradiance and quasinormal modes. 
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